Crowd synchrony and quorum sensing in delay-coupled lasers 
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Crowd synchrony and quorum sensing arise when a large number of dynamical elements communi- 
cate with each other via a common information pool. Previous evidence in different fields, including 
chemistry, biology and civil engineering, has shown that this type of coupling leads to synchroniza- 
tion, when coupling is instantaneous and the number of coupled elements is large enough. Here we 
consider a situation in which the transmission of information between the system components and 
the coupling pool is not instantaneous. To that end, we model a system of semiconductor lasers 
optically coupled to a central laser with a delay. Our results show that, even though the lasers are 
non-identical due to their distinct optical frequencies, zero-lag synchronization arises. By changing 
a system parameter, we can switch between two different types of synchronization transition. The 
dependence of the transition with respect to the delay-coupling parameters is studied. 

PACS numbers: 05.45.Xt, 42.65.Sf, 42.55.Px 



Many situations in nature involve systems of multiple 
dynamical elements that interact with each other through 
a common medium. Examples include pendulum clocks 
mounted on the same wooden beam [1], cellular pop- 
ulations communicating via small molecules that freely 
diffuse among the cells [2, 3], and longitudinal modes 
of a laser connected through saturation of the common 
amplifying medium [4, 5]. In most cases, the coupling 
surrounding leads to a synchronous behavior between all 
the coupled elements, but no general framework for such 
crowd synchrony has been developed yet. 

One of the features that determines whether a coupling 
medium is able to elicit synchronization is the number of 
elements that are connected through it. One situation 
where this fact has been established is that of the Mil- 
lennium Bridge in London. Two days after its opening 
in 2000, the pedestrian bridge had to be closed because 
of excessive wobbling. Subsequent investigations [6] re- 
vealed that the pedestrians, who initially walked with 
different frequencies and phases, fell into a synchronous 
pacing when the number of people on the bridge was 
large enough. That effect was modeled and understood 
in terms of the synchronization of simple oscillators [7, 8]. 
In a biological context, collective glycolytic oscillations in 
yeast cells have also been seen to arise as a result of large 
enough cell densities [9, 10]. In that case, and in con- 
trast with the Millenium bridge, oscillations of individ- 
ual cells only occur above a critical number of them: the 
coupling induces both the oscillations and the synchro- 
nization among them. A similar situation was recently 
reported on a chemical oscillator system, formed by cat- 
alytic particles suspended in a catalyst-free Belousov- 
Zhabotinsky medium. This system also shows a transi- 



tion to synchronized oscillations above a certain density 
of catalytic particles [11]. Furthermore, depending on 
the coupling strength, the synchronization appears either 
progressively or suddenly. The question of how generic 
this behavior is remains open. 

Studies of crowd synchrony to date have considered 
the coupling with the medium to be instantaneous [7- 
11]. This naturally results in synchronous behavior with 
zero time lag between any pair of elements in the system. 
However, in many situations the transmission of the cou- 
pling signal takes an amount of time that is non- negligible 
with respect to the characteristic time scales of the sys- 
tem components. This is the case e.g. in systems of 
technological importance such as optically coupled semi- 
conductor lasers. When these devices are separated dis- 
tances on the order of centimeters, they are subject to 
coupling delays on the order of the characteristic time 
scales of solitary lasers (shorter than nanoseconds). In 
recent years much effort has been devoted to understand 
the synchronization of semiconductor lasers. This is im- 
portant for technological reasons, i.e., to achieve large 
output powers, but also for increasing our knowledge of 
how generic dynamical systems synchronize. Semicon- 
ductor lasers are low-cost, versatile, and many of the 
commercial lasers are well-described theoretically. They 
also show a wide variety of nonlinear dynamical behav- 
ior, both as single elements with external influences and 
as part of laser networks; examples include low- frequency 
fluctuations [12, 13], chaos [14], excitability [15-17], and 
nontrivial synchronization phenomena [18, 19]. However, 
most studies of coupled lasers so far have considered a 
small number of elements. Thus, how to achieve syn- 
chronization for a large number of coupled non-identical 
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lasers is still an open question. 

In this Letter we show that a collection of M semi- 
conductor lasers coupled through a central laser exhibits 
zero-lag crowd synchronization. This setup is a gener- 
alization of the case of two identical oscillators coupled 
through a third central element [20-25]. Isochronal syn- 
chronization is relevant in both technological [26] and 
biological [27] contexts. Here the central laser operates 
in a passive regime (below threshold), and plays the role 
of a coupling medium analogous to the bridge structure 
in pedestrian synchronization [7], and to the reaction 
medium in chemical synchronization [11]. Our results 
show that the general properties of both the crowd syn- 
chrony and the quorum-sensing transition are readily re- 
produced with this setup. 

The equations describing the slow envelope of the com- 
plex electric field E and the carrier density TV for the M 
lasers are [23] 
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+ e \Ej,h\ ), and the subscripts H 
and j denote the central (hub) laser and jth outer (star) 
laser, respectively. The field and carrier decay rates are 7 
and 7 e , respectively, a is the linewidth enhancement fac- 
tor, e is the gain saturation, ujq is the optical frequency 
and ujh and Uj are the detunings of the hub and the star 
lasers with respect to the reference frequency uoq. The 
coupling is characterized by its strength n and delay r, 
both of which are assumed equal for all lasers, fij = ji 
and fin are the pump currents of the star and hub lasers, 
respectively. Finally, £j(t) and are uncorrelated 

complex Gaussian white noises, with D being the noise 
strength. The model was integrated with the stochastic 
Heun algorithm with a time step of 0.8 ps, using parame- 
ter values typical for semiconductor lasers: 7 = 300 ns -1 , 
7 e = 1 ns -1 , a = 3. D = 10 -5 ns -1 , ujq = 2ttc/X (where 
c is the speed of light and A = 654 nm). ujh — without 
loss of generality, and ujj is chosen from a Gaussian dis- 
tribution with zero mean and standard deviation a. In 
what follows we neglect nonlinear gain saturation (e = 0), 
since it does not affect the results obtained (not shown). 
Unless otherwise stated a = 207T rad/ns, corresponding 
to AA =0.014 nm, and r = 5 ns. 

Figure 1 shows the stationary emitted intensity for 
varying number of star lasers. For small M [Fig. 1(a)] 
the lasers oscillate independently. By increasing M, syn- 
chronized emission at near zero-lag occurs for lasers with 
similar frequencies, forming clusters with similar dynam- 
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FIG. 1: (Color online) Synchronization of star-coupled semi- 
conductor lasers for increasing number of elements, M. The 
color coding shows the intensity for each star laser as a func- 
tion of time. In the vertical axis the lasers are sorted by their 
solitary frequency, ujj, with number 1 corresponding to the 
most negative detuning. The black dashed line shows the 
laser for which cj = 0. (a) M — 10, (b) M 25, (c) M = 30 
and (d) M = 75. The right column shows the frequency 
ujj/27T of the lasers (dots), in relation with the normailized 
cumulative Gaussian distribution (solid line) . The pump cur- 
rents are \i — 1.02 and \ih = 0.4, and the coupling strength 
k = 30 ns -1 . 



ics as shown in Fig. 1(b). The number of synchronized 
lasers in those clusters grows as M increases [Fig. 1(c)], 
with an emission characterized by short pulses of irregu- 
lar amplitudes with a repetition period around 2r. Those 
characteristics become more evident for large M, where 
almost all the lasers emit synchronously at zero lag [Fig. 
1(d)], with emission pulses taking place simultaneously 
in most of the lasers. This situation is reached provided 
the pump current of the hub laser is set below the solitary 
lasing threshold, i.e. \±h < l^th — 1 [28]. 

In order to quantify the level of zero-lag synchroniza- 
tion, we calculated the total coherent intensity of the star 
2 
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Figure 2 compares the dynam- 
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ics of this quantity with that of the intensity \Eh\ 2 of the 
hub for increasing number of lasers. For increasing values 
of M, Figs. 2(a-d) shows the emergence of large intensity 
pulses in the total coherent intensity, corresponding to 
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FIG. 2: (Color online) (a-d) Time trace for the intensity of the hub laser, \Eh\ 2 (top trace, blue), and for the coherent intensity, 
/ (bottom trace, black), (e-h) Cross-correlation function between \Eh\ 2 and /. The number of lasers is M — 25 (a,e), M — 30 
(b,f), M = 75 (c,g) and M — 100 (d,h). The parameters are the same as in Fig. 1. 



strongly synchronized activity in Figs. 2(b-d). The hub 
laser reproduces these dynamics after a time r. This is re- 
flected in a large peak at time — r in the cross-correlation 
function between \Eh\ 2 and /, shown in Figs. 2 (e-h). 
Thus, the hub laser lags behind the star lasers in the 
synchronized state. 

To investigate the transition to the synchronized state, 
we use as order parameter the time-averaged total coher- 
ent intensity of the star lasers (/), where (•) is the aver- 
age over a time window of length T = 4 /is, calculated in 
the stationary state. In the absence of synchronization 
(/) grows linearly with M, while when synchronization 
arises this linear dependence is lost. Figure 3(a) shows 
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FIG. 3: (Color online) Ratio between the averaged coherent 
intensity (/) and the number of star lasers M, as a function 
of M itself and for different coupling strengths: k — 10 ns _1 
(red), k — 20 ns _1 (blue) and k — 30 ns _1 (black), (a) \i — 
1.02, [ih = 0.4. (b) [i = 0.99, [ih = 0.4. (c) fi = 0.7, [ih = 0.4. 
Each point is averaged over 10 to 40 different initial conditions 
and detuning frequencies. The arrows mark errorbars out of 
the axis limits. 

the average coherent intensity as a function of the number 
of star lasers for different coupling strengths and pump 
currents. When the star lasers are pumped above the 
solitary threshold and for small coupling, (I) /M is ap- 
proximately constant, corresponding to the case of the 
absence of synchronization. For moderate values of ft, on 
the other hand, the system becomes gradually synchro- 



nized as M increases. The transition to synchronization 
occurs for a critical number of coupled lasers M c , which 
we quantify as the number of lasers for which the growth 
rate of (/) with M changes abruptly. For even larger ft 
the critical number of lasers needed for synchronization 
decreases. 

The qualitative character of the synchronization tran- 
sition can be changed by tuning the pump current /i of 
the star lasers below the laser threshold. When \i is well 
below threshold [Fig. 3(c)], the transition to synchroniza- 
tion is very sharp, in contrast with Fig. 3(a) above, pro- 
vided coupling is large enough. Note that in this case 
both the star and hub lasers are pumped below their 
solitary threshold, but are effectively above threshold due 
to coupling, and they turn on due to their spontaneous 
emission. The transition takes place when the star lasers 
are pumped at their solitary threshold [Fig. 3(b)], which 
shows the beginning of a sharp transition for intermediate 
ft (blue line online) when the star lasers are pumped only 
slightly below threshold. We also note that this behavior 
requires that the hub laser be pumped below threshold, 
i.e. fiH < l^th = 1, otherwise the critical behavior is lost. 

One of the characteristic features of crowd synchro- 
nization is the scaling of the critical number of ele- 
ments with the degree of diversity in the population and 
with the coupling coefficient [7]. In our case diversity is 
caused by the different frequencies ujj of the lasers. Fig- 
ures 4(a,b) show the dependence of the critical number of 
lasers on the width a of the frequency distribution and 
on the coupling strength ft. The results show that M c 
increases linearly with a, while the dependence with ft 
follows a power law with negative exponent, as occurs 
in [7]. As expected, the larger a the more different the 
lasers, and more lasers are required to generate the syn- 
chronized state. A broad frequency distribution leads to 
a reduction in the size of the clusters of lasers with simi- 
lar ujj showed in Fig. 1(b), and a corresponding decrease 
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in the coherent intensity. On the other hand, the larger 
the coupling strength the smaller the minimum number 
of lasers required to synchronize the system [Fig. 4(b)]. 
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FIG. 4: (Color online) Critical number of lasers, M c , as a 
function of different system parameters, for pump currents 
ji = 1.02, \±h = 0.4. (a) M c as a function of the width of the 
frequency distribution, a. The dashed line shows a linear fit 
for ct/2tt > 6 GHz (M c oc 1.2a, r 2 = 0.982). Other parame- 
ters are k = 30 ns -1 and r = 5 ns. (b) Doubly logarithmic 
plot of M c as a function of the coupling strength k. The 
dashed line shows a power-law fit of the data (M c oc 1/ft 2 ' 2 , 
r 2 = 0.987). Other parameters are a — 20tt rad/ns and r = 
5 ns. (c) M c as a function of the time delay, r. Other param- 
eters are k — 30 ns -1 and a = 20tt rad/ns. The simulations 
are averaged over 20 stochastic realizations of the initial con- 
ditions and frequency distribution. 

We have also considered the effect of the coupling delay 
r on the transition to the synchronized state. As shown 
in Fig. 4(c), for short delays (compared with the charac- 
teristic time scales of the laser) M c exhibits a sharp reso- 
nance at a r corresponding to the cavity decay time, but 
for longer delays M c is reduced and is almost indepen- 
dent of r. When the coupling delays are not identical [29], 
results similar to those of Fig. 3 are found, but for larger 
coupling strengths. In that case the synchronized dy- 
namics may be characterized by periodic fluctuations of 
small amplitude, or even steady state emission. 

In conclusion, we have shown numerically that a sys- 
tem of non-identical semiconductor lasers coupled to a 
common hub laser with time delay can be synchronized 
with zero lag. The transition to the synchronization oc- 
curs above a certain critical number M c of coupled lasers, 
provided the pump current of the hub laser is smaller 
than the solitary pump current threshold fi t h- The type 
of synchronization transition can be controlled via the 
pump current of the star lasers: a gradual (second-order- 
like) transition is observed for star lasers with \i > /j, t h, 
and an abrupt (first-order-like) transition arises for \i < 
/ith' A similar behavior has been exhibited by a chem- 
ical quorum sensing system [11]. The critical number 
of lasers increases linearly with the width of frequency 
distribution, and depends on the coupling strength via 
a power-law with negative exponent, in agreement with 
the crowd synchronization transition reported in the Mil- 
lenium bridge [7]. On the other hand, the coupling delay 
reduces the critical number of lasers while it has no in- 



fluence on it for large enough time delays, even though 
the delay is evident through the lag time with which the 
hub laser is synchronized with the star lasers (which are 
synchronized isochronously to one another). 
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